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Abstract 

This work is devoted to the study of discrete ambiguities. For parametrized potentials, they 
arise when the parameters are fitted to a finite number of phase-shifts. It generates phase equiva- 
lent potentials. Such equivalence was suggested to be due to the modulo n uncertainty inherent to 
phase determinations. We show that a different class of phase-equivalent potentials exists. To this 
aim, use is made of piecewise constant potentials, the intervals of which are defined by the zeros 
of their regular solutions of the Schrodinger equation. We give a classification of the ambiguities 
in terms of indices which include the difference between exact phase modulo n and the numbering 
of the wave function zeros. 



pacs: 24.10.-i ; 03.65.-w ; 03.65.Nk 



1 Introduction 



Attempts to determine the potential from scattering data have a long history, and conditions to obtain 
a unique answer are well known (see for instance the textbooks by Newton [T], Chadan and Sabatier 
[2J). Approaches to the three-dimensional inverse scattering problem can be classified in two categories 
[3 [5] . In the first case, known as fixed-£ problem, the potential can be constructed from the phase- 
shifts 5(1, fc), if they are known for all momenta k g (0, oo), and from the discrete spectrum (bound 
state energies and the corresponding normalization constants). Note that the potential is assumed to 
satisfy adequate integrability conditions [2]. In the second case, the so-called fixed- E problem, Loeffel 
[3] has obtained theorems ensuring a unique potential from the knowledge of the phase-shifts 5(i, k) 
for all (non-discrete) non negative values of A = £ + 1/2. If the data set reduces to discrete values 
of A = £ + 1/2 for non-negative integer £, the Carlson's theorem [J| predicts a unique potential V(r), 
provided it belongs to a suitable class [H [3] . 

The present work is dealing with the second aspect of the problem. It is relevant to the case 
where a chosen parametric form is used to fit a differential cross-section at a fixed energy. Here, 
the uniqueness is defined in the sense of a best fit to a differential cross section, for instance. The 
argument does not apply to phase shifts, which can be extracted from the scattering amplitude, since 
a phase is defined mod nir. Such undeterminations are called phase- or discrete ambiguities. They 
have been noticed long ago by Drisko, Satchler and Bassel [5]. 
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As a concrete example, let us quote the discrete ambiguities described by Goldberg et al [BJ. 
By using potentials of Woods-Saxon shape, these authors have extracted four a — 208 Pb potentials, 
characterized by depth differences of 50 McV in the real potential part, giving similar ^-square per 
degree of freedom. The x' s square are related to the phase-shifts via the fit to the differential cross 
sections. The Jeffreys, Wentzel, Kramers and Brillouin (JWKB) phase-shifts have been calculated 
for the real parts of the four potentials. They were found to differ by a factor mr, n integer and 
independent of £ for low enough angular momentum (£ < 44). 

Attempts to study these ambiguities are scarce. Sabatier [7J and Cuer [8J have studied their origin 
by means of the JWKB approximation. A physical interpretation has been given by Leeb and Schmid 
[9 J , in which the occurrence of discrete ambiguities is linked to the existence of partly Pauli forbidden 
states. 

Note that the discrete ambiguities are still of interest. Let us mention the recent papers on the 
subject by Brandan [TO] and by Amos and Bennett [TTJ. These works have been done in the context 
of heavy and light nuclear ion scattering, respectively. 

The purpose of the present work is a further study of the discrete ambiguities. In order to discuss 
this problem in detail, and stress the origin of the modulo n shift, we shall consider piecewise-constant 
real potentials. Such potentials have been used recently by Ramm and Gutman |12j . These authors 
have shown that different piecewise-constant positive real potentials can lead to almost the same 
fixed-energy phase-shifts. As we shall see, this is a convenient starting point for the present study. 
Here we enlarge the study and consider also negative potentials. Moreover, we found it particularly 
useful to restrict the class of potentials to the subset of piecewise-constant potentials the intervals of 
which are defined by the zeros of the regular solution of the Schrodinger equation. 

The present work shows that piecewise constant potentials are a basis for a class of phase equivalent 
potentials, with a phase ambiguity of mr. Within each n, these potentials can be ordered according 
to the positions of the zeros of the regular solutions. The value of n has a minimum, which can be 
calculated by means of the JWKB approximation. 

Dealing with the exact solution of the Schrodinger equation, our work goes beyond, for a certain 
class of potentials, the earlier attempts by Sabatier [7\ and Cuer [8j, who have resorted to the JWKB 
approximation. The drawback is that our method becomes very tedious as the number of partial 
waves increases. 

In what concerns numerical applications, the energy will be put equal to unity, without loss of 
generality. Indeed, use is made of the following scaling property: if 5(£,k = 1) is the phase for the 
potential V(r) at the energy E = 1, it is also the phase for the potential k 2 V(kr) at the energy 
E = k 2 . 

The paper is organized as follows. In section 2 the basic formalism is recalled. In section 3 our 
results are presented. The results are discussed in the JWKB approximation in section 4 and our 
conclusions are presented in section 5. 



2 Formalism 

In this paper we investigate to what extent piecewise constant potentials allow us to reproduce a 
finite set of fixed-energy phase-shifts. Considering spherically symmetric potentials of finite range, we 
start from the reduced radial Schrodinger equation. For a positive energy E = k 2 and an angular 



2 



momentum I (H = 2m = 1), it reads 

MK r)" + (k 2 - v(r) - ^J^) r)=0. (1) 

Here, the reduced potential is denned by 

w(r ) = 2 -^ V (r) , (2) 

in terms of the potential V(r), and tpg(k,r) is the regular solution of Eq.([T]) defined by the Cauchy 
condition \im r ^o ■ipe(k,r)r~ e ~ 1 = 1. 

Our construction of phase-equivalent potentials is based upon the following property. Suppose the 
function ipg(k,r) to vanish for r = r(£,k). It is unnecessary to know the potential for r < r(£,k) 
to determine the phase shift S(£,k). This holds for (£,k) fixed. Thus, it is sufficient to consider 
the Schrodinger equation on the interval [r(£,k),oo[. Actually, the derivative ip' e with respect to r 
is undetermined at r = r(£,k), which has no consequence since the phase-shift is determined only 
by the ratio ipg/ipg. Furthermore, for a fixed energy and a fixed £, there is a countable number of 
zeros [13] • The functions tpe(k,r) and ip' e (k,r) cannot vanish at the same value of r for 2£ + 1 > 
(except for r = 0) and I > as shown in [14] . Consequently, ipi(h,r) has only simple zeros (except 
at r = 0). We recall that thanks to the fact that the zeros are simple, they satisfy the strict ordering 
< h(i,k) < r 2 {£,k) < ... < f n (£,k) <...<R. 

The nth zeros of the regular solution satisfy a monotony property, namely, for fixed n, the function 

£^f n (£,k) (3) 
increases with t as has shown by Sturm in the 1830s [TS]- More precisely we have 

i-Q- r Mk,f n (£,k)) 

which is positive definite. As stated before, tp^ and ip 1 * cannot vanish simultaneously for 2£ + 1 > , 
except for r — and £ > [14] , so that the denominator is never equal to zero. 

To find potentials reproducing N + 1 fixed energy E — k 2 phase-shifts, S(£j, k),£ = 0, 1, . . . N, we 
apply the following procedure : 

1. We construct a piecewise constant potential with a compact support: v(r) — for r > R, R 
being fixed. 

2. f n (£,k) is the n-th zero of the regular solution ipi{k,r) of the Schrodinger equation (JTJ) for the 
£-wave aside from the trivial zero at the origin. We denote by rg a generic non trivial zero of 
tpg(k,r). For the zero potential f n (£,k) is given in terms of the zeros of the spherical Bessel 
function ji(z) labelled X£ >n . We have kr n {£, k) = Xi >n . 

3. The potential v of Ea.([2"]) is assumed to be constant on each interval delimited by the partition 
< ri < r 2 < r 3 , . . . < r N < R. We set 

v(r) — vo, K 2 = k 2 — vq ; < r < r% 
v(r) = i>i, K\ = k 2 — v\ ; r\ < r < r 2 

v(r) = vn, K 2 N = k 2 - vn ; r N < r < R (5) 
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4. The regular ue(kr) and irregular wi(kr) solutions of Eq. ([T]) for v = are denoted, respectively, 

wi{x) = -J-^-Y e+1/2 (x) 

in terms of the Bessel functions J^,!^ of order v, given in [16] , We have ut(x) = xje(x) where 
ji is the spherical Bessel function of order t, 

5. Let ip£ be the regular solution of the Schrodinger equation ([lj. We denote by A(r), B(r), C(r), D(r), 
the value of the ratio 

^(k,r) W 
for I = 0, 1, 2, 3, ... respectively. Here, the prime denotes the first derivative with respect to r. 



As we shall see, these requirements give rise to a class of phase-equivalent potentials. 

To show phase-ambiguities, it is worth to recall that the exact phase can be calculated by using the 
variable phase method of Calogero [T7]. Most of the time it allows the removal of the nir ambiguity 
burdening the direct solution. With this method, the phase-shift is reached by solving a first order 
differential equation 

^-8(£, k,r) = -^p- ( Ui (kr) cos{5{£, k, r)) + w e {kr) sm(S(£, fc, r))f , (7) 
or k 

with 8(i, k, 0) = as boundary condition. The phase-shift is given by the limit 8(1, k) = liim.^.,^ 5(1, k, r). 
In the present work, use is made of the variable phase method to determine the exact phases. In 
other words, k being fixed, we solve Eq.([7]) for each value of r starting from r — 0, combined with 
8(1, k, 0) = 0. Actually, lim^^ 5(1, k, r) = 5(£, k, R) = 8(1, k). 

A variant of the Calogero's equation but for tan(S(£, k,r)) has been derived in [18]. It can be 
obtained from Eq.© noticing that 

|tan(^, fc ,,)) = cos( ^ fc;r))2 ls(£,k,r). (8) 
Setting T(£, k, r) = taa(5(£, k, r)), it yields 

^T(£, k,r) = -^ (ui(kr) + w e (kr)T(S(£, k, r))) 2 , (9) 

which is solved with the boundary condition T(l, k, 0) = 0. At zero energy, the latter formalism allows 
the calculation of both the scattering length and the effective range directly from the potential. In 
our case, we have chosen to work with Eq.© to avoid the numerical difficulties owing to poles in 
T(£,k,r). 



3 Examples 

3.1 A single phase-shift 

Although the case of a single known phase-shift is trivial, it illustrates features of the discrete ambi- 
guities. 
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Assume, for example, the s-wave phase-shift 8{£ = 0, k) = So to be known. The procedure amounts 
to find a constant vq in [0, R\. Setting 

tan(fc f + ' o) = A(R) , (10) 
k 

the equation to be solved is, depending upon whether K 2 is positive or negative 

ta,n(K R) 
K~ 

tanh^i?) 
K* 



= A(R) , K = JK$ 



A(R), K* =J-K*. (11) 



For > 0, discrete ambiguities come from the fact there is a countably infinite set of solutions 
to Eq. lfTTj) . Whatever the right-hand-side of (jlll) . one solution of (fTTj) exists in the interval 

(2n - 1)| < \/k 2 - vo R < (2n + 1)| (12) 

for every n > 1. In [0,7r/2[ the solution exists only for suitable values of the r.h.s. of the equation. 
The difference between the solutions for Vq on two adjacent intervals is of the order of the square of 
the length of the interval, namely n 2 . At the potential level, it introduces a difference of AV ~ 2mR 2 ' 
As a consequence, the larger is the mass m, the smaller are the differences between phase equivalent 
potentials, in agreement what is observed phenomenologically [19] ■ It means that while searching 
optical potential parameters, the chance to find phase equivalent potentials on a finite parameter 
space is larger for heavy than for light systems. Indeed discrete ambiguities are often observed for 
nucleus- nucleus scattering with a heavy target, for instance for a — 208 Pb scattering (not so for 
a — 58 Ni), whereas few attempts have been made for nucleon-nucleus scattering. Note that for proton 
scattering on Cr isotopes at 10 MeV, discrete ambiguities have been reported by Andrews et al. [20J . 
The difference in the potential depth is about 100 per cent. At high energies I is weak as compared 
to fc 2 , and we expect \/k 2 — voR to fall in the same interval as kR, namely 

(2n fe -l)| < k R < (2n fc + l)| . (13) 

This is due to the fact that \Jk 2 — vq~R ~ kR — voR/(2k). As a consequence, the solution cannot be 
found in another interval, and the discrete ambiguities disappear. This is confirmed by the analysis 
of scattering data at high energies. Finally, to finish with this example, we have verified that the 
different solutions of Eq. (fT2|) lead to s-wave phase shifts which differ by multiples of w. 



3.2 2 phase-shifts £ = 0, 1 



The determination of the piecewise constant potential starts from the largest zero rjv < R fitting the 
solution vn{Kn) on the largest £ — N phase-shift considered. Finding vn-i(Kn-i) on the previous 
interval [rN-i,rN[ is more complicated in general, as it depends not only on rjv-i but also from 
rjv, Kj\f . The complexity increases with the number of zeros considered. Two exceptions exist, namely 
the £ = 0, 1 and £ = 0, 1, 2 cases. These exceptions are called simply soluble, and their solutions are 
presented below. 

The regular solution for £ = 1 is denoted u\(Kor) for r < r\. Since r\ is a zero of ipi, we have 
K Q rx = xi tTll , and 

(Vni) wi(xi, ni ) = sin(a; l! „ 1 )/x 1:ni - cos(x^ ni ) = . (14) 
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On the interval [ri, R[, up to a constant factor, it reads : 

ipi(k,r) = ui(Kir)wi(Kin) - ux{Kiri)wx{Kxr) . (15) 

Finally, for r > R the regular solution is proportional to 

tpi(k, r) — ui(fcr) cos(5i) + w\(kr) sin(^i) . (16) 

Equating the inverse of the logarithmic derivative of solutions (I15I16P at r = R, we obtain 

tantfiTiCR-n)) Af 

K[ = V Kl 

tanh(K^(R-ri)) Af 

K* ~ V 

Here we have 



Kl = J-K\ . (17) 



Af = —B{R)R -R 2 + B(R)n + Rn - B(R)K 2 R 2 n 

V = —B(R) — R + B(R)K 2 R 2 - B(R)K 2 Rri — K\R 2 r\ . (18) 

We recall that in these expressions, B(R) = (ipi/tfj[)(k, R), with ipi given by Eq. tfTcj]) . at r > R. 

Turning to the I = wave, the regular solution in [0, ri[ is proportional to sin(ifor). Therefore, 
(tpo/tfj' )(r = n) = tan(_Ko?"i)/-Koj which is equal to r\ when Eq. (jT^J) is taken into account. At the 
other end of the interval, [r±, R[, the function ipo/ip' i s given by 



ip (k, R) _ Kxtx cos(isTi(J? - n)) + wx(Ki{R - n)) 
%(k,R) ~ Ki{cos(Ki(R - n)) - K in sm(Ki(R - njj) ' 

It is identified with A(R) — taa(kR + 6o)/k. We then obtain 

taxi(Ki(R - n)) A(R) - n 



Ki 



(19) 



Kl = J-Kl . (20) 



Ki ~ 1 + A(R) Kl n 

tanhiKfjR-n)) A(R) - n 

K* 1 + A(R) K\ ri 

Equating equations (IT71) . (US]) andf |2"0)) , we are left with 

a _ A(g)g(g) + A(g)ij - g(g)ij - i? 2 

1 A(R)B(R)R 2 ' 1 j 

Provided that A(R)B(R) ^ or equivalently that i? is a zero of neither ipo nor The equation (|2"TT) 
also reads 

o o Af 

K\ = k 2 + — 
A/" = fc 3 .Rsin(£ -£r) 

2? = sin(fc-R + <5 ) (sin(fei2 + 5i) - kRcos(kR + 5i)) . (22) 

It fixes the value of K 2 , and consequently the constant vi, on the interval [ri, i?[ independently of the 
value of ri. For this very reason, the model is quoted as simply soluble. With this value of K\ 1 Eq. 
([20| allows us to determine n . In this case, a first source of ambiguities come from the choice of the 
value of n among the zeros of ipi smaller than R. The next step consists in determining vq. To this 
aim, use is made of the fact that KqTi — x\^ ni is a zero of ipi Vni. The choice of which zero xi^ ni of 
the spherical Bessel function ji is another source of ambiguity. 
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x x ,i = 4.493409 


x h2 = 7.725252 




n 


K 


n 


K a 




-7 


0.467351105 


-6 


0.803489067 


9.61463257 


-6 


0.524561926 


-5 


0.901848241 


8.5660223 


-5 


0.597764867 


-4 


1.02770172 


7.51701782 


-4 


0.694775396 


-3 


1.19448618 


6.467427 


-3 


0.829515446 


-2 


1.42613676 


5.41690836 


-2 


1.02946743 


-1 


1.76990236 


4.36479003 


-1 


1.35771439 





2.33423792 


3.30953931 





1.99999964 


1 


3.4384809 


2.24670498 


1 


3.89058288 


2 


6.68884865 


1.15494498 



Table 1: For the first two zeros of ji(x), the value of Kq is listed as function of the shift (nw) with 
respect to the exact phase of the reference potential Kq = 2, r\ = 2.24670498, K\ = 3,R= 10. 

Numerical application 

To fix ideas, the following example is treated. A reference potential is defined by setting vq — 
—3, r € [0, r%[, vi — —8, r G [ri,R[. The energy is fixed at E = k 2 = 1. Thus, we have Kq — 2 and 
K\ = 3. The first zero of the spherical Bessel function j% is chosen, namely x\ t \ = 4.493409. It fixes 
the value of T\ at 2.2467045. The two phase-shifts So and Si are then calculated by the variable phase 
method. 

With these pseudo-data, we look for phase-equivalent potentials of the same range R — 10, i.e. 
potentials giving the same £ = 0, 1-phase-shifts modulo 7T. In fact the phase-equivalent potentials are 
determined from the values of A(R), B(R), R = 10 calculated from the reference potential. We recall 
that K\ is recovered from Eq. ([21]) . and is therefore identical to the starting value. In the Table 1 are 
shown the different values of r% < R = 10 obtained from Eqs. ([20|) . Only the first two zeros of the 
spherical Bessel function ji(z) = u\(z)/z namely 2:1,1 = 4.493409, xi. 2 = 7.725252, are displayed in 
Table 1. They concern the ambiguities relating to Kg. 

The ambiguities proceed from two sources: the different values of ri and the choice of the zero 
xi t m- As a consequence, two different potentials may give the same shift with respect to the exact 
phase. At this stage, the ambiguities can be characterized by two numbers: 

1) ri\ : numbering the zero of the regular solution for the p-wave (£ = 1) 

2) n : the difference between the exact phases of the phase-equivalent and the starting potentials 
divided by tt. In the present example this number happens to be the same for both I = and £ = 1 
phases. 

The values of K are listed and ordered as function of the mi shift with respect to the exact 
solution. For n — and X\ t i = 4.493409, the parameters of the reference potential are recovered 
with a 10~ 7 accuracy. Due to the finite range R, Table 1 exhausts the number of phase equivalent 
potentials for the first two zeros of j\. Other cases exist, corresponding to the higher zeros of j%. Here 
n\ = 1 for all A'o of the second column and n\ = 2 for the A"o's of column 4. 

It is interesting to quote results for transparent potentials. They are characterized by Sq = S\ = 0. 
It implies automatically K\ = k = \f~E~ = 1 and n takes only the two values 7.72525183 or 4.49340947 
smaller than the range parameter R — 10. A freedom still exists in the choice of Xi tni which determines 
Kq. The results are displayed in Table 2. 
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n 


m 


K 


ri 


-1 


l 


0.581652127 


7.72525183 





l 


0.999999931 


4.49340947 





2 


1.00000004 


7.72525183 


1 


2 


1.71924064 


4.49340947 


1 


3 


1.41149086 


7.72525183 


2 


3 


2.42669234 


4.49340947 


2 


4 


1.82080713 


7.72525183 


3 


4 


3.13040547 


4.49340947 


3 


5 


2.22915123 


7.72525183 


4 


5 


3.83244722 


4.49340947 


4 


6 


2.63697599 


7.72525183 


5 


6 


4.53359608 


4.49340947 


5 


7 


3.0444899 


7.72525183 


6 


7 


5.23421053 


4.49340947 



Table 2: Transparent potentials R — 10. For two values of n = #1,1,0:1,2, Kq is listed as a function 
of n defined in the text, together with m, the order of the zero of ji(x). 

3.3 3 phase-shifts £ = 0, 1, 2 

Consider the case where three phase-shifts 60,61,62, corresponding to £ = 0, 1 and 2, are known. On 
the interval [r 2 ,i?[, the I = 2 regular solution reads, up to a constant multiplicative factor, 

i\>i{k,r) = U2{K 2 r)w2{K 2 r2) - ^(i^r^w^i^r) . (23) 

It vanishes at r = r 2 as imposed. For r > R, it is proportional to 

i>2(k,r) — u 2 (kr) cos(5 2 ) + w 2 (kr) sm(6 2 ) ■ (24) 

With C(R) = (ip2 / ^'2) (k,R) ^ the continuity conditions at r — R require 



tan(K 2 (R - r 2 j) 
K 2 

ta,nh(K^(R- r 2 )) 
K* 

Here we have 

TV = -18C(R)R-9R 2 + 3C{R)KiR 3 + l8C{R)r2+9Rr 2 - 9C(R)Ki R 2 r 2 

- 3K 2 R 3 r 2 + 6C(R)K^Rrl + ZKlR 2 rl - C{R)K^R 3 rl 

V = -18C(R) -9R + 9C{R)KlR 2 + SK^R 3 - l8C{R)KlRr 2 - 9K\R 2 r 2 
+ 3C(R)K 2 l R 3 r2 + 6C(R)K 2 r 2 + 3K 2 Rr 2 - 3C(R)K%R 2 r 2 

- KlR 3 r 2 . (26) 



±, K 2 
V 

* K * 

V ' 2 



(25) 



Other constraints come from the fact that r 2 is a zero of ip2- First, for < r < r\, the function 
ip2 is proportional to u 2 - At r = n, we have (u2/u' 2 )(K ri) — — ri/2, since for x = K ri we have 
sin(a;)/a; — cos(ir) = u\(x) — 0. Expanding -02 as a linear combination of U2{Kir) and w 2 (ifir) (in 
the case where K\ > 0,Ki ~ y/K 2 ) on [ri,r 2 [, taking into account that (7/> 2 /-0 2 )(fc, r) = —r\/2 at 
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r = r\ and zero at r — r 2 , we are left with: 

tan(ifi(r 2 - n)) _ TV 

K x V 

tanh(JTf(ra-n)) _ TV 



if* 2? 

TV -3n + 3r 2 + iffrirl 



(27) 



2? 3 + 3Kfnr 2 - Kfr. 

For ^ = 1. the regular solution is zero at r = r±, and reads at r = r 2 , up to a non zero multiplicative 
factor, (when iff > 0, K x = y/Kf) 

ipi(k,r 2 ) = Ui(Kxr 2 )wi(Kxri) - wi{Kxr 2 )u\{Kirx), (28) 

and its derivative is 

ip[{k,r 2 ) = K^u'^Ktr^iKm) - w'^Kxr^miKm)) . (29) 

It implies that 



r 2 (ri + T+r 2 (-1 + K\ r\ T)) ^\%r 2 ) 



-Ti+n-T- K\ ri r 2 r + K 2 r\ (n + T) V4(&, ra) ' 

Here we have defined 

tan(ATi(r 2 - ri)) A" " 

r = ^ *i = V*i 

tanh(if 1 *(r 2 -n)) 



(30) 



r = ^ Kt = yJ-iq. (31) 

Equating T with TV/2? of Eq.(27J| and using Eq.p?|). we obtain [ipi/i/>i)(k, r 2 ) = B(r 2 ) = r 2 /2. 

The ^ = 1 regular solution ipi on [?"2,2?[ is a linear combination of u\(K 2 r) and u^i^O (when 
2ff — 0: Ki = >/ Ki ) ■ Using the constraint that (^i/^)(fc, r) is equal to r 2 /2 at r = r 2 and to 23(2?) 
at r = 2? we obtain: 



tan(^ 2 (^-r 2 )) TV 

iV2 



2f 2 2? 
tanh(2f 2 *(2?- r 2 )) _ TV 



K* = J-Ki (32) 



TV = W(R)R + 32? 2 - 3B{R) ri - 32?r 2 

+ 3B(R)K%R 2 r 2 - B{R)KlRr\ - K\R 2 r\ 

V = 3B(R) + 32? - ZB{R)KlR 2 + 3B(R)K%Rr 2 

+ 3K\R 2 r 2 - B(R)K\r\ - K\Rr\ 

+ B{R)K\R 2 r\ . (33) 



Identifying ([25j [26]) and (|32I33|) we have finally 



2 _ 4i?(2?)C(2?) + 2B(R)R - 2C(R)R - R 2 
K2 ~ B(R)C{R)R 2 (M) 
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k z 



or, equivalently, provided that R is not a zero of tpi or tp 2 

a 

V 

N = k 5 R 3 sin(<5i - S 2 ) 
V = {sm.(kR + 8 1 )-kRcjos{kR + 5 1 )) 

x ((3 - k 2 R 2 ) sin(fci? + 8 2 ) - 3kRcos(kR + 8 2 )) 



(35) 



Here again the value of K 2 and of v 2 on [r 2 , R[ does not depend on r 2 , and the model is still simply 
soluble. Once K\ is known, we have to determine r 2 through Eqs. ()32|33j) . Then Eq. (f2T1l expresses 
K 2 in terms of the ratio (ipe/ip' e )(k, r 2 ) for I = 0, 1, namely the quantities A{r 2 ),B(r 2 ). We recall that 
(tfji/^p' 1 )(k,r 2 ) = r 2 /2 — B(r 2 ), as shown before. By definition A(r 2 ) — (iljo/ip' )(k,r 2 ). Thus, on the 
interval [r2,-R[, the matricial equation for K 2 > 



cos(K 2 (R-r 2 )) sMK 2 (R-r 2 )) 
-K 2 sin(K 2 (R-r 2 )) cos(K 2 (R - r 2 )) 



ipo(k,r 2 ) 
%{k,r 2 ) 



(36) 



allows the determination of A(r 2 ) (for K\ < the extension of the matricial equation is straightfor- 
ward) . Indeed, the solution of this matricial equation yields 



A(r 2 ) = 



A(R) - T 



with 



T 



A(R)K 2 T + 1 

t&n(K 2 (R - r 2 )) 
K 2 

ta,nh(K%(R-r 2 )) 



K-2 



K* 



K* 



(37) 



(38) 



To determine K 2 on the interval [ri , r 2 [ use is made of Eq. ([21]) by setting R — r 2 , namely 

, A(r 2 )B(r 2 ) + A(r 2 )r 2 - B(r 2 )r 2 - r 2 
K i = 



where B(r 2 ) = r 2 /2. 
Let us define 



A(r 2 )B(r 2 y 2 



F Q = UA{R)C{R) + 27 A(R)R - 5AC(R)R - 27R 2 

-27 A(R)C(R)R 2 - 9A{R)R 3 + 9C(R)R 3 + mA{R)C{R)r 2 2 
+lSA{R)Rrl - mC{R)Rrl - l8R 2 r\ + 18C(i?)r| + 9Rr\ 



(39) 



(40) 
(41) 



F, 



A8A{R)C{R)R 2 r 2 - 6A(R)R 3 r^ + 6C(R)R 3 r 2 
-18A(R)C(R)Rrl + 9A(i?)i? 2 ^ - 9C{R)R 2 rl - 3R 3 r 

F 3 = -ZA{R)C{R)R 3 r 3 2 



(42) 
(43) 



Gi = -18A(R)C(R)r 2 - 9A{R)Rrl + 18C(R)Rr 2 2 + 9R 2 r\ 
-18C(R)rl - 9Rr 3 2 



(44) 



G 2 = 



§A(R)C{R)R 2 rl + ZA{R)R 3 rl - 3C{R)R 3 r% - \8A{R)C{R)Rrl 
-9A{R)R 2 r 3 2 + 9C{R)R 2 rl + 3R 3 r 3 2 + 6A(R)C(R)4 
+3A(R)Rr 2 - 6C(R)R4 ~ 3i?M 



(45) 
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G 3 = +SA(R)C(R)R 3 r 3 - 3A{R)C(R)R 2 r* - A{R)R 3 r 2 i + C(R)R 3 4 . 



(46) 



We get the following equation 



F + F x Kl + F 2 Kl + F 3 Kl + G X K\ + G 2 K\kI + G z K\K\ = . 



(47) 



It gives K\ in terms of A(R) , B (R) , C (R) and K\,r 2 . Then r\ is determined (0 < r\ < r 2 ) through 
Eq. (f2"7| . The value of Kq is such that KqT\ = X\ >ni , whatever n\, forcing r\ to be a zero of 

For scattering data corresponding to a potential in the considered class, namely constant on inter- 
vals determined from zeros of the regular solution, the phase-equivalent potentials all correspond to a 
permutation of the zeros of the regular solution with the constraint that (W = 1, . . . , N — 1) r^, n< < 

n+i, ne+1 and r N ^ nN < R. 

The transparent potentials are zero on [ri,i?[. The values r\ and r 2 are respectively zeros of the 
spherical Bessel functions j\ and j%. The potentials are non-zero only on [0,ri[. The value of Ko, 
such that Kori = x\ differs from k = y/E if r\ ^ x\/k. This is what leads to the discrete ambiguities. 

Numerical Application 

We construct a reference potential, which is constant on each interval [0, r\[, [ri,r2[, [ra,i2[. Here, 
we choose vq = —3,Kq = 2,vi = — 8, Ki = 3,v 2 = —4:,K 2 = y/Z, together with r\ = 2.24670598 and 
7~2 = 2.6958027, which reduce ipi and ip 2 to zero, respectively. The phase-shifts Sj,j = 0, 1, 2 of this 
potential are calculated at E = 1 by means of A(R), B(R), C(R). 

Phase equivalent potentials are found by fitting piecewise constant potentials to these phase-shifts 
or more explicitly A(R), B(R), C{R). First, the values of K\ and v 2 are equal to the starting values 
v 2 = -4, if 2 = \/5 according to Eq. ([34]) or Eq. ((35j) . Then, r 2 is determined from Eqs. (|32I33|) 
with the constraint that r 2 < R. The value of K\ is obtained from Eq. (f39|) . and r\ is determined 
(0 < r\ < r 2 ) through Eq. (|27|) . The value of Kq is such that Kqt\ — ri being a zero of tpi. 

Once the phase equivalent potential is determined, the corresponding Schrodinger equation is 
solved for each partial wave. It yields the respective positions of the zeros of the regular solution. The 
resolution of Eq. ([7]) yields a phase-shift from which one subtracts the phase-shift of the reference 
potential. 

A sample of solutions for R = 5.5 and R — 10.0 are displayed in Tables 3 and 4, respectively. The 
results underline a classification according to three indices : 

- n, the difference between the phase calculated directly from Eq. ([7]) and the exact value divided 
by 7r. It is found to be the same for the three waves £ = 0,1,2. 

- Tlx, numbering the zero of regular solution of the £ = 1 wave, 

- n 2 , numbering the zero of the regular solution of the £ = 2 wave. 

The reference potential is recovered for the {n, n\, n 2 } = {0, 1, 1} case. Moreover, in each subclass 
of fixed n, the solution can be ordered according to {rii, n 2 } in such a way that n 2 > n\. 
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Table 3: Phase-equivalent potential to the reference potential of line 2. The range parameter is 
R = 5.5. Kj — k 2 — Vj, k = 1. The classification is made according to indices n, m, ni (see text) 
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Table 4: Same as Table 3 for a range parameter R — 10.0. The reference potential appears on line 5. 
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3.4 4 phase-shifts I = 0, 1, 2, 3 

We begin the procedure by determining the value of K3 as a function of r^. On the interval [r3,i?[, 
the regular solution ^3, which vanishes at r = 7-3, is a linear combination of 113 (K 3 r) and w^{K^r) : 

ip3(k,r) = u 3 (K 3 r)w 3 (K 3 r 3 ) - u 3 (K 3 r 3 )w 3 (K 3 r) . (48) 

According to the definition ([S]), the ratio (ip 3 /ip' 3 )(k, R) = D(R) at r = R. It yields 

tan(K 3 (R-r 3 )) AT 3 



K 3 V 3 
tanh(J£S(ii-r3)) _ A/" 3 



= : i iff = »/-J« , (49) 

K 3 V 3 3 V 3 ' v ; 

where Af 3 ,T> 3 are given in Appendix A. We first check if the potential is simply soluble. Looking at 
Eqs. (|21l34p . we expect 

2 _ 9C(R)D(R) + 3C(R)R - 3D(R)R - R 2 

Ks ~ C(R)D(R)R 2 (50) 
or, equivalently, provided that R is not a zero of i\) 2 or ip 3 



K$ = k 



2 



V 



Af = k 7 R 5 sin(S 2 -63) 

V = ((3-k 2 R 2 )sm(kR + 5 2 )-3kRcos(kR + 8 2 )) 

x (3(5 - 2k 2 R 2 ) sm(kR + S 3 ) + kR(k 2 R 2 - 15) cos(fci? + 83)) ■ (51) 

The equation (|50p is verified if and only if (ip 2 /ip' 2 )(r = r 3) = C( r 3) = r 3/3- This can be proved in 
the following way. On the interval [r 3 , R{, we express ip 2 as 

Mk,r) = a(K 3 )u 2 (K 3 r) + (3(K 3 )w 2 (K 3 r) . (52) 

The ratio a(Ks) / (i(K3) is determined in terms of C(r 3 ) and C(R), which yields 

taa{K 3 {R-r 3 )) = ^3 
K 3 V>3 



K3 = JK § 



tanh(if 3 *(E-r 3 )) N'3 _ , „ 2 
K* = — K 3=\/~ K 3> 



K* 3 =J-KI, (53) 



'3 u 3 

where M' 3 and T>' 3 are given in Appendix A. 

By equating N3/V3 = Af' 3 /T>' 3 , and taking iff from ([50]). we find indeed C(r 3 ) = r 3 /3. 

It remains to actually calculate C(r 3 ). To do so, we first note that D(r{) = 5r\/{x\ n — 15). Here, 
use is made of the fact that K^ri = xi. ni is a zero of ipx and thus sin(A'o r i) = K^r\ cos(ito?"i)- 

On the interval [7*1 , r 2 [, ip3 is written as 

i>3{k, r) = aiKx^Ktr) + p{K 1 )w3{K l r) . (54) 

The ratio a{Ki)/f3{Ki) is given as a function of £)(ri), thus as a function of ri,xx tni . This ratio 
being known, D(r 2 ) can be calculated. 

Taking into account Eq. (|27p . D(r 2 ) is rational fraction in terms of x% tni , K 2 , n, r 2 . 
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Finally, on [r 2 , r 3 [, ip 3 is given by 

i> 3 (k, r) = a(K 2 )u 3 (K 2 r) + P(K 2 )w 3 (K 2 r) . (55) 

Similarly to the preceding step, the ratio a(K 2 )/ f3(K 2 ) is a function of D(r 2 ). Inserting its value 
in (|55p . and by using Eqs. (|25[) and (f2"6")l at i? = r 3 , we obtain a compact polynomial expression 
involving X\ tni , K 2 , K 2 ,rx,r 2 ,r 3 ,C{r 3 ). By remembering that Z?(r 3 ) = and setting C(r 3 ) = r 3 /3, 
the condition for a simple soluble model reads 

(9 + 3K 2 r 2 + #M)(225 + 45K 2 r 2 + %K\r\ + K & 2 rl){K 2 r 2 - x\ ni ) = . (56) 

Setting K\r\ = 3x, K 2 r 2 = y and ifi n = z, the polynomial in a; has no real root, the polynomial 
in y has a single real root, namely K\r\ = —5.39254. The polynomial in z has root K\T\ = x\^ ni 
(K 2 > 0) implying Kq = K%. Consequently, C(r 3 ) = r 3 /3 occurs only in exceptional situations, and 
the model is not simply soluble. 

Since the model is not simply soluble we let r 3 run from to R. Use is made of Eq. (|49]) and 
the first equation of Appendix A to determine K 3 . The value K 3 {R — r 3 ) is determined in intervals 
]m 3 n, (m3 + l)7r[, where m 3 (m 3 > 0) is fixed at each run. We have ensured that a single solution exists 
in each interval, with K 3 being a continuous function of r 3 and ip 3 (r 3 ) — 0. Then A(r 3 ), B(r 3 ),C(r 3 ) 
(see Eq. ©) for £ — 0,1,2, respectively, are determined from A(R) , B (R) , C (R) by using the fact 
that on [r 3 , R[ we have 

ipi{k,r) = a e u e (K 3 r) + f3 e wi(K 3 r) . (57) 

By comparison with the previous case of 3 phase shifts, for which K\,K 2 are independent on r\,r 2 
respectively, many more solutions have to be considered. 

We first look at transparent potentials, for which the exact phase is nir for all waves £ = 0,1, 2, 3. 
The results are displayed in Tables 5 and 6 for R = 12 and R = 15, respectively. 

We denote by n the class of phase-equivalent potential, for which the exact phase is nir for all 
waves £ — 0,1, 2, 3. The present study emphasizes these classes. Each class is characterized by the set 
of numbers {m, n 2 , n 3 } where ng is n^-th zero of the regular solution of the Schrodinger equation for 
the £-wave. The configurations can be ordered according to {ni,n 2 ,n 3 } with n\ < n 2 < n 3 < n max . 
The number of configurations and n max for fixed n increase with R. 

Obviously, except for the lowest exact phase, different sets of indices n\ < n 2 < n 3 exist, with 
different Kj,Vj giving the same phase nir. Consequently, the exact phase is not a sufficient criterion 
to remove the phase-ambiguities, as suggested in ref.[5_. The size of the set associated to the index n 
increases with the value of n. 

Obviously, for each n, the set of minimal size corresponds to {ni, 712,^3} = {1, 1, 1}; [n max = 1). 
For the lowest n, this set is the only solution. Note that in Table 5 and 6 for n = 0, all the solutions 
reduce to the same zero potential. 



As a more realistic example, the reference potential is chosen as the starting potential with pa- 
rameters E = k 2 — 1, R = 10, 



Vfc 2 


- vq 


= K 


= 2 


< r < n = 2.24668994 




- Vi 


= K X 


= 3 


n < r < r 2 = 2.699578766 


Vk 2 


- v 2 


= K 2 


= v/5 


r 2 <r <r 3 = 7.38655698 


Vk 2 


- v 3 


= K 3 


= y/2 


r 3 < r < 10 . 



1 Note that here we may encounter transparent potentials for which the fourth exact phases are n(£)ir, with n(£) not 
the same for all the waves. 
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n 
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K x 
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1.000 
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1.000 
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1.000 
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4.493 


1.000 


9.095 


1.000 


10.417 


1.000 


12.000 





2 


2 


2 


1.000 


7.725 


1.000 


9.095 


1.000 


10.417 


1.000 


12.000 


1 


1 


1 


1 


1.141 


3.936 


1.249 


4.968 


1.483 


5.716 


1.310 


12.000 


1 


1 


1 


2 


1.171 


3.838 


1.251 


4.864 


1.398 


8.081 


1.240 


12.000 


1 


1 


1 
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1.174 


3.827 


1.237 


4.863 


1.340 


10.681 


1.191 


12.000 


1 


1 


2 


2 


1.155 


3.892 


1.264 


7.530 


1.538 


8.316 


1.313 


12.000 


1 


1 


2 


3 


1.161 


3.870 


1.240 


7.579 


1.423 


10.728 


1.233 


12.000 


1 


1 


3 


3 


1.160 


3.872 


1.272 


10.021 


1.599 


10.812 


1.318 


12.000 


1 


2 


2 


2 


1.119 


6.901 


1.355 


7.933 


1.765 


8.616 


1.421 


12.000 


1 


2 


2 


3 


1.135 


6.805 


1.292 


7.881 


1.537 


10.784 


1.288 


12.000 


1 


2 


3 


3 


1.126 


6.863 


1.367 


10.231 


1.853 


10.907 


1.430 


12.000 


1 


3 


3 


3 


1.079 


10.102 


1.878 


10.889 


3.562 


11.232 


2.029 


12.000 



Table 5: Transparent potentials such that the first four phase-shift are all equal to nir. K 2 = 
k 2 — Vj\ k = 1. The range of the potential R = 12. The index ne is such that re is the ne-th zero of 
the regular solution in the £-wave. 



The functions A(R) , B (R) , C (R) , D(R) are calculated, and the procedure developed in subsection 3.4 
is used to determine phase-equivalent potentials. 

The classification of the solutions is performed very much in the same way as for transparent 
potentials. Here n is the difference between the exact phases of the reference and the phase-equivalent 
potentials divided by ir. A priori, n depends on £. As an example, the following phase-equivalent 
potential vj = k 2 — K 2 

K ,n = 3.12729631,1.43683507 
K 1 ,r 2 = 1.39601562,4.82494432 
K 2 ,r 3 = 2.54346784,4.97266696 

K 3 ,R = 0.820438746,10 (59) 

is such that the first 3 values of n{£ = 0) — n(£ — 1) = n(£ = 2) — -2, while n(t = 3) = -3. In this 
case, the 3 indices ne are not ordered according to n± < n 2 < n 3 . Actually, n\ = l,n 2 = 2,n 3 = 1. 



Let us concentrate on potentials for which n is independent of £. Results are displayed in Table 7. 
The same conclusions apply to this ensemble as to the case of transparent potentials. We note that 
one of the n = solutions reproduces the reference potential. 



The configuration n = 1 for non transparent potentials is reported in Appendix B. 



4 Generalization by means of the JWKB approximation 



The method developed in the preceding sections is not suitable for an extension to a larger number of 
phase shifts, as it would require a tremendous numerical effort. In order to work beyond 5 phase-shifts 
and obtain some general results, use can be made of the JWKB approximation. We recall that in this 
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Table 6: Same as Table 5 but for R = 15 
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Table 7: Same as Tables 5 and 6 but for a non zero starting potential exhibited between the double 
lines. 
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case the phase-shift is given by: 

8{£,k) = lim f / K(£,k,r') dr' - f K(£,k,r')dr') (60) 

r ^ co \Jrt(g,k) Jft(Lk) J 

(remember that E = k 2 ). For the type of potential considered in this work, we have 



K(£,k,r) = \jk 2 -v{r)- li * l + l) k(£,k,r) = y ' k 2 - ^ ? + ^ . (61) 

In (|60p. rt(£,k) is "the" largest turning point relative to the function if considered, while rt(£,k) — 

^W+T)/k. 

The JWKB phase 5(£,k) has been calculated for the transparent potentials of Table 5, and com- 
pared to the exact phase. In case of the null potential, the JWKB estimate is obviously exact. For 
non zero potentials, the relative differences 

_ s(e,k)-S(e,k) 

D *- ] — wj) — 1 ( } 

are most of the time better than 1 %, but as large as 2 % in a few instances. 

Furthermore, the JWKB method can be shown to predict the position of the zeros of the regular 
solution. To this aim, the regular solution is written 



^(fc,r) ocsin ( / \/k*-v{r')- e ( £ ^ dr' ) . (63) 

rt(£,k) 

Clearly, for f n (£, k) to be the nth zero of the regular solution, the following condition has to be 
satisfied : 

ffnim I £(£+1) 



k 2 — v(r') — — dr = mr . (64) 

rt(£,k) 



Actually, the quantity to be considered is 



h = -( r (e ' k) J k 2 - v(r>) ~ e -^^ dr') + eg . (65) 

I" \Jrt(£,k) V r' 2 I 



The levelling value eg is given by 



c< = 1 - - J^P*', (66) 

where r±{£) is the first zero of the spherical Bessel je(z). The quantity eg does not depend on the 
value of A;. It is designed in a way that here Ig is an integer and yields the exact result 1 for the first 
zero in the absence of potential. 

Calculations have been performed for the transparent potentials of Table 5. The ratio between the 
exact values and the JWKB estimates show an excellent agreement. The deviation with respect to 
unity is most of the time better than 1 %, with few exceptions of the order of a few %. In this case, 
Ig is close to an integer and represents the position of the nth zero of the regular solution for the £ 
wave. 

Comforted by the good quantitative agreement of the above calculations, the JWKB approximation 
can be used to ask questions on a general basis. Note that TV > 3 is required, otherwise the model is 
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simple soluble and the JWKB approximation is not necessary. For instance, it is possible to predict 
the lowest value of n for a set of N + 1 phase-shifts, and a fixed potential range R. According to the 
results displayed in the previous Tables, it occurs for the configuration in which all the rijsr = 1. Thus, 
it is sufficient to consider the largest value of N. 

The data are the number of phase shifts N + 1 the energy E = k 2 and the potential range R. 
Recalling the previous conventions, fi(N,k) denotes the first zero of the regular solution for the 
iV-wave, supposed to be strictly lower than R. 

We first consider transparent potentials. By definition, we have 

Q1 + Q2- Q3 

n = , (67) 

with 

rfi{N,k) j N(N + 1) 

Qi = / \/fc 2 - v(r') - { / > dr' , (68) 

Jrt(N,k) V r u 



Q2 = l R Jk 2 -v{r>)-- {N + l) dr' , (69) 



and 



f 

JJn(n+i 



0,-1 M-amn*. ( 70 ) 



y/N(N + l)/k 

The first contribution yields 

Ql = 7T - cntt ■ (71) 
However, c/v is defined by Eq. (|66[) . which yields 



rri(N)/k 

Qi= , \/fc 2 - ""V - dr' . (72) 

Jy/N(N+l)/k 



-V) k v r 



The potential is constant on the interval [fi(N, k), R[ so that the two next integrals are analytical and 
read 

Q, = JkI& - N(N + 1) - ^N(N + 1) aictan J + L) (73) 

- v^w ij - m + 1) - smrr^un ^gj«Effl ± 9 

and 



g 3 = y/k*R*-N(N + l) v^V+lJarctan ^ ^^^±lj . (74) 

It is clear that for the configurations having all n^r = f , the value of n is ordered by the value of 
fi(N, k). The lowest configuration corresponds to fi(N, k) being the closest to R. The absolute limit 
is given, provided kR high enough and R ^ f\(N, k), by 

n m in > — — , (75) 

because Q2 is strictly positive. In a compact form taking into account the equation (|72[) we obtain 



* J ri {N)/k V r 
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Actually, n m - m is integer and we have to take 

l ri (N)/k 



> 1 



k2 _ N{N+1) 



(77) 



Here, [x] denotes the integer part of x, i.e. the integer m x such that m x < x < m, + 1. In the present 
case we use the property that any integer m which satisfies m > x satisfies m > [x] + 1. 

If we take as an example the upper case of Tables 5 and 6, the numerical estimate of Eq. ([75|) yields 
"-min > —1-47631 and n m i n > —2.39885 for R = 12 and R — 15, respectively. The integer values 
satisfy n m - m > — 1 and n m - m > —2, respectively, in agreement with the values of Tables 5 and 6 for 
R = 12 and R = 15. 



On the other hand, by looking at the Tables 5,6 and 7, the maximum value of n 3 depends on n. 
It obeys the following recurrence relation 

i3,max(n + 1) = rc. 3 , max (n) + 1 ■ (78) 



For the transparent potentials, n3. max (0) is given by the number of zeros of the regular spherical 
function j 3 (kr) below kR (we exclude the exceptional situation where a zero of j 3 is equal to kR). For 
instance let us consider R = 12, k = 1. There are two zeros of j 3 below 12, namely 6.988 and 10.417. 
They can be found in the Table 5 concerning the configuration n = (zero potential). This implies 
that 7i 3 ,max(0) = 2 so that rc. 3 , max (-l) = n 3)max (0) - 1 = -1 (see Table 5). 

A similar example is found from Table 6. Three zeros of j 3 are found below kR = 15 : 6.988, 10.417 
and 13.698. It yields fi3,max(0) — 3, which implies that 7i3 )m ax(— 2) — 1. We recover the lowest 
configuration of the Table 6. 

Yet a further example is provided us by R = 21, k = 1. Five zeros of j'3 exist below 21 : 
6.988,10.417,13.698,16.924 and 20.122. It gives rc 3>max (0) = 5 and n 3 , max (-4) = 1. This config- 
uration has been confirmed numerically. It corresponds to the potential (vj = k 2 — K 2 = 1 — K? ): 

K a ,ri = 0.436410526,10.2962892 
Kx,r 2 = 0.34393,13.836594 
K 2 ,r 3 = 0.321277,20.0371883 

K 3 ,R = 0.791895048,21. (79) 
The exact phases of this potential are — 4ir for I = — 3. 



Taking the same case, R = 21, k = 1, but for 5 phases instead of 4, 4 zeros of are situated below 
21, namely 8.183,11.705,15.040 and 18.301. As expected, n 4 , max (0) = 4 and n 4 , roax (-3) = 1. The 
lowest configuration is given by the following potential (vj = k 2 — K 2 = 1 — K 2 ): 





n 


= 0.536260007,8.37916127 


K x 


r 2 


= 0.218548173,11.03235 


K 2 


ra 


= 0.360733527, 13.5849 


K 3 


'4 


= 0.793613806,14.7336 


K 4 ,R 


= 0.266900301,21. 



(80) 

Its exact phases are — Sir for £ = — 4 



Beyond 5 phases, finding numerical solutions becomes too tedious to be interesting. 
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Taking into account the results for 4 and 5 phases (£ — 0, 1, 2, 3 and I = 0, 1, 2, 3, 4), we propose 
the following conjecture : 



Consider models which are not simply soluble, i.e. number of phases N + 1 strictly greater than 
3, and configurations such that n is the same for all partial waves, the minimal value of n is given by 
"min = ^ n {jN) + 1) where n{jjq) denotes the number of zeros of the regular spherical Bessel function 
jN strictly smaller than kR, kR not being a zero of j'jv and high enough to ensure n(jj\r) =/= 0. 

As stated in the preceding section, our attention has been focused on cases such that the phase 
difference is mr, independently of I. However, phase equivalent potentials with n depending on £ can 
be constructed. Two examples are quoted below for transparent potentials. 

The first example corresponds to 4 phase shifts, namely Sq = 5\ = S 2 = 2tt and S 3 = tt. The phase 
equivalent potential reads 

K ,n = 6.80886924,0.659934689 
Ki,r 2 = 1.2151,12.8731444 
K 2 ,r 3 = 1.34685,13.8434964 

K 3 ,R = 1.15320976,15. (81) 
The values of ng are and n\ — 1, n 2 = 5, n 3 = 4. 

The second case includes 5 phase shifts, namely 5q = Si = 62 = tt and ^3 = ^4 = 0. The phase 
equivalent potential is given by 

K ,n = 7.43592354,0.604283917 
K u r 2 = 1.08584,5.38409355 
K 2 ,r 3 = 1.13832,12.2708388 
K 3 , ri = 1.07046,13.398 

K 4 ,R = 0.814872,21. (82) 
The values of rig are n\ = 1, n 2 = 2, n 3 = m = 3. 



A similar bound is valid for non transparent potentials. We remind the reader that in this case 
n7r is the difference between the exact phase of the equivalent potential and the exact phase of the 
reference potential. For non transparent potential the above inequality becomes 



u min > - (tt-cntv- [ Jk 2 -v(°)(r>)- JV(JV /2 +1) rfr / | . (83) 

\ Jrt(N,k) V r' 2 J 

Here, is the starting reduced potential and rt(N, k) the corresponding largest turning point. On 
the basis of the results obtained for transparent potentials, we have 




2 - N ^^ dr'+d(N,k)j . (84) 

The quantity S(N, k) is either the exact phase of the starting potential calculated from Eq. (JTJ) or the 
JWKB approximated phase Eq. (|60p. 

For N — 3 and the starting potential of Table 7 we obtain n min > 1 + [-4.08774] = -4 and 
"-min > 1 + [—4.01251] = —4 according to the fact that the phase is exact or approximated. It has to 
be compared to the value n m ; n = —3 of Table 7. 
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Let us make a brief comment about the maximal number of nodes N of the wave functions (except 
the trivial zero at the origin ) inside the potential. For a given £, N is the number of zeros ri of 
the regular solution of the Schrodinger equation below R. We have checked, from examples of Tables 
5-7 and Appendix B, that for all potentials in the "class" n (such that (5(£, k) = nir,£ = 0, 1, 2.3) N 
depends only on £, k, R, n. Moreover N(£, k, R, n) satisfies N(£, k, R, n+ 1) = N(£, k, R, n) + 1. In this 
sense, n is correlated to N as expected from the results of Sabatier [7|, who argued that additional 
multiples of ir were associated with additional nodes "inside" the potential. However the Sabatier's 
method, based on the Abel transform, cannot be directly applied here. Indeed the latter method 
requires the function r H> r 2 (k 2 — v(r)) to be continously differentiable [8]. 



5 Conclusions 



The present work is devoted to the study of discrete ambiguities. They arise when a specific parametrized 
expression is used as a potential, its parameters being fitted to a finite number of phase-shifts, or more 
directly to a scattering amplitude. To this aim, use is made of piecewise constant potentials, the in- 
tervals of which are defined by the zeros of the regular solution of the Schrodinger equation. These 
potentials generate a class of phase equivalent potentials with a phase ambiguity of mr, n being an 
integer. The reference potential belongs obviously to this class. 

A few examples have been investigated and solved exactly for 1-5 phase-shifts. The number of 
discrete ambiguities depends on the range of the potential. For a given value of n, the phase equivalent 
potentials can be ordered according to the position of the zeros of the regular solutions, except for 
the zero potential. 

Note that the present study considers both positive and negative potentials. We have verified that 
in the case of attractive potentials having at least one bound state, the measurement of the ground 
state energy is a sufficient criterion to fix the potential and get a unique answer. 

For a number of phase-shifts larger than 5, the numerical effort becomes rapidly intractable for an 
exact solution. In this case, use can be made of the JWKB approximation. This last allows to derive 
a lower bound for the minimal value n m j n . 

Advantage has been taken at several places of the transparent potentials. They are defined as 
potentials having all their phase-shifts equal to nir. Particular attention has been given to the trans- 
parent potentials for which n is independent of £. Cases can occur, however, such that the exact 
phases are not the same for all waves. Then n depends on £. 

Moreover, for transparent potentials, a conjecture has been proposed to determine n m \ n . 

Finally, we show that removing the phase-ambiguity as suggested by Drisko.Satchler and Bassel 
is not sufficient to determine the potential uniquely. Potentials reproducing the exact phases can 
differ by the numbering of the zeros of their regular solutions. Actually, the difference between the 
two approaches lies in the following. If a scattering amplitude is reproduced with a Woods-Saxon 
potential, for instance, the number of parameters is limited, and the large number of phase-shift to 
be fitted ensure the quasi unique determination of the parameters. Thus, phase ambiguities can be 
observed, and possibly removed. On the other hand, piecewise constant potentials generate a much 
large functional space. Consequently, ambiguities arise also from the choice of the intervals. This 
behavior has been emphasised by Lombard and Wilkin |21| . in analysing high energy differential cross 
sections via the Glauber model. In this case, there is no phase ambiguity but the potential is defined 
as the statistical average over a large ensemble of different piecewise constant potentials fitting the 
data. 
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Appendix A 



In this Appendix we report the expression of J\f 3 ,T> 3 of Eg. p5|) and J\f' 3 ,V 3 of Eg. (IS"5)) 

A/3 = 675D{R)R + 225R 2 - 9QD{R)KlR 3 - IbK^R 4 - 675D(R)r 3 

- 225Rr 3 + 315D(R)K 2 R 2 r 3 + 90K 2 R 3 r 3 - 15£>(i?)A^i? 4 r 3 

- 27QD{R)KlRrl - 90K%R 2 r 2 + 36D(R)KlR 3 r 2 
+ 6K 4 R 4 rl + A5D{R)K 2 r 3 + l5K 2 Rr 3 

- 21D(R)KlR 2 rl - 6AT^i? 3 r| + D(R)K 3 R 4 r 3 

V 3 = Q7hD{R) + 225R - 315D(R)K 2 R 2 - 90K 2 R 3 + lhD{R)K 4 R 4 

+ 675D(R)K 3 Rr 3 + 22hK\R 2 r 3 - 90D(R)K 3 i R 3 r 3 - 15AT| R 4 r 3 

- 27QD{R)Klr 2 3 - 9QKlRr 2 3 + 126 D(R)Kl,R 2 r 2 3 
+ mK 4 R 3 rl - 6D(R)KlR 4 rl - AbD{R)K 4 Rrl 

- 15K 3 i R 2 4 + &D(R)KlR 3 r 3 3 + KlR 4 r 3 3 



A/"' 3 = 9R 2 n - 9Rr 2 3 + SK^R 3 ^ - 3K 3 2 R 2 4 + l8Rr 3 C{R) - 3K^R 3 r 3 C{R) 

- l8r 2 C(R) + 9K 2 R 2 r 2 C(R) - 6K 2 Rr 3 C(R) + K 4 R 3 r 3 C(R) + \8R 2 C{r 3 ) 

- l8Rr 3 C(r 3 ) + 6K 2 R 3 r 3 C(r 3 ) - 9K 2 R 2 r 2 3 C(r 3 ) + 3K 2 Rr 3 C(r 3 ) 

- K 4 R 3 r 3 C{r 3 ) + 36RC(R)C{r 3 ) - 6K 2 R 3 C(R)C(r 3 ) - 36r 3 C(R)C(r 3 ) 
+ \8K 2 R 2 r 3 C{R)C{r 3 ) - 18K 2 Rr 2 C(R)C{r 3 ) + 3K 4 R 3 r 2 C(R)C(r 3 ) 

+ 6K 2 r 3 C(R)C(r 3 ) - 3K 4 R 2 r 3 C(R)C(r 3 ) 

V' 3 = 9Rr 3 - 3KlR 3 r 3 + 9KlR 2 rl ~ SK^Rr^ + K^R 3 ^ + 18r 3 C(R) 

- 9K 2 R 2 r 3 C(R) + l8K 2 Rr 2 3 C(R) - ZK 4 R 3 rlC{R) - QKjr 3 C(R) 

+ 3K 4 R 2 r 3 C(R) + 18RC(r 3 ) - QKlR 3 C{r 3 ) + \8KlR 2 r 3 C{r 3 ) - 9K 2 Rr 2 3 C(r 3 ) 

+ 3K 4 R 3 r 2 C(r 3 ) - 3K 4 R 2 r 3 C(r 3 ) + 36C(R)C(r 3 ) - 18K 2 R 2 C(R)C(r 3 ) 

+ 36K 2 Rr 3 C(R)C(r 3 ) - 6K 4 R 3 r 3 C(R)C(r 3 ) - l8K 2 r 2 C(R)C(r 3 ) 

+ 9K 4 R 2 rlC{R)C{r 3 ) - 6KlRr 3 C(R)C(r 3 ) + KlR 3 r 3 C{R)C{r 3 ) 
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Appendix B 



n 


m 


n 2 


n 3 


K 


n 


K x 




K 2 


r3 


K 3 


R 


1 


l 


1 


1 


1.737 


2.587 


3.446 


2.990 


2.883 


3.118 


2.456 


10.000 


1 


l 


1 


2 


0.918 


4.896 


3.900 


5.275 


1.855 


7.097 


4.680 


10.000 


1 


l 


1 


3 


1.932 


2.325 


3.475 


2.722 


2.673 


5.356 


2.209 


10.000 


1 


l 


1 


4 


1.996 


2.252 


3.437 


2.650 


2.593 


6.627 


2.075 


10.000 


1 


l 


1 


5 


2.057 


2.184 


3.361 


2.589 


2.499 


8.024 


1.895 


10.000 


1 


l 


2 


2 


1.943 


2.312 


3.476 


3.635 


2.737 


3.805 


2.197 


10.000 


1 


l 


2 


3 


2.106 


2.134 


3.515 


3.442 


2.596 


4.903 


2.019 


10.000 


1 


l 


2 


4 


2.235 


2.011 


3.505 


3.309 


2.473 


6.180 


1.832 


10.000 


1 


l 


2 


5 


2.360 


1.904 


3.456 


3.233 


2.339 


7.675 


1.595 


10.000 


1 


l 


3 


3 


2.171 


2.070 


3.492 


4.302 


2.549 


4.518 


1.886 


10.000 


1 


l 


3 


4 


2.388 


1.882 


3.539 


4.085 


2.373 


5.716 


1.639 


10.000 


1 


l 


3 


5 


2.592 


1.734 


3.539 


3.938 


2.187 


7.225 


1.332 


10.000 


1 


l 


4 


4 


2.395 


1.877 


3.475 


5.035 


2.309 


5.309 


1.506 


10.000 


1 


l 


4 


5 


2.648 


1.697 


3.520 


4.818 


2.068 


6.729 


1.139 


10.000 


1 


l 


5 


5 


2.514 


1.787 


3.367 


5.991 


2.009 


6.339 


1.032 


10.000 


1 


2 


2 


2 


2.032 


3.802 


4.212 


4.148 


2.772 


4.323 


2.384 


10.000 


1 


2 


2 


3 


1.324 


5.834 


5.330 


6.117 


2.084 


7.739 


4.610 


10.000 


1 


2 


2 


4 


2.228 


3.467 


4.100 


3.821 


2.573 


6.558 


2.032 


10.000 


1 


2 


2 


5 


2.327 


3.320 


3.974 


3.682 


2.441 


7.909 


1.783 


10.000 


1 


2 


3 


3 


2.208 


3.498 


4.069 


4.634 


2.589 


4.852 


2.000 


10.000 


1 


2 


3 


4 


2.378 


3.248 


4.027 


4.395 


2.428 


5.986 


1.745 


10.000 


1 


2 


3 


5 


2.544 


3.037 


3.923 


4.189 


2.238 


7.393 


1.418 


10.000 


1 


2 


4 


4 


2.389 


3.234 


3.912 


5.225 


2.336 


5.500 


1.565 


10.000 


1 


2 


4 


5 


2.615 


2.954 


3.867 


4.967 


2.094 


6.856 


1.182 


10.000 


1 


2 


5 


5 


2.483 


3.112 


3.691 


6.078 


2.019 


6.429 


1.054 


10.000 


1 


3 


3 


3 


2.346 


4.647 


4.825 


4.956 


2.639 


5.171 


2.126 


10.000 


1 


3 


3 


4 


2.480 


4.398 


4.691 


4.714 


2.492 


6.258 


1.869 


10.000 


1 


3 


3 


5 


2.628 


4.150 


4.490 


4.479 


2.303 


7.581 


1.530 


10.000 


1 


3 


4 


4 


2.488 


4.382 


4.479 


5.418 


2.369 


5.692 


1.631 


10.000 


1 


3 


4 


5 


2.689 


4.056 


4.326 


5.127 


2.125 


6.989 


1.231 


10.000 


1 


3 


5 


5 


2.551 


4.274 


4.111 


6.170 


2.033 


6.524 


1.079 


10.000 


1 


4 


4 


4 


2.657 


5.293 


5.134 


5.586 


2.401 


5.858 


1.693 


10.000 


1 


4 


4 


5 


2.833 


4.966 


4.864 


5.273 


2.155 


7.109 


1.280 


10.000 


1 


4 


5 


5 


2.685 


5.240 


4.592 


6.252 


2.046 


6.608 


1.103 


10.000 


1 


5 


5 


5 


2.858 


6.026 


5.112 


6.322 


2.059 


6.678 


1.124 


10.000 



Table 8: This Table is an extension of Table 7 (see text). 
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